Incoherent matter-wave solitons 
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The dynamics of matter-wave solitons in Bose- Einstein condensates (BEC) is considerably affected 
by the presence of a surrounding thermal cloud and by condensate depletion during its evolution. 
We analyze these aspects of BEC soliton dynamics, using time-dependent Hartree-Fock-Bogoliubov 
(TDHFB) theory. The condensate is initially prepared within a harmonic trap at finite tempera- 
ture, and solitonic behavior is studied by subsequently propagating the TDHFB equations without 
confinement. Numerical results demonstrate the collapse of the BEC via collisional emission of atom 
pairs into the thermal cloud, resulting in splitting of the initial density into two solitonic structures 
with opposite momentum. Each one of these solitary matter waves is a mixture of condensed and 
noncondensed particles, constituting an analog of optical random-phase solitons. 



The physics of quantum-degenerate, interacting bose 
gases closely resembles the behavior of light in non- 
linear media. The dynamics of Bose-Einstein conden- 
sate (BEC) at zero-temperature is well-described by the 
Gross-Pitaevskii (GP) mean-field theory, given by the 
nonlinear Schrodingcr equation (NLSE) for the conden- 
sate order parameter. The same equation describes the 
evolution of coherent light in nonlinear Kerr media. This 
analogy has opened the way for the field of nonlinear 
atom optics |l|, |2| with striking demonstrations of famil- 
iar nonlinear optics phenomena such as four wave mix- 
ing , superradiant Rayleigh scattering Q , and matter- 
wave amplification |£J. |(|. carried out with matter- waves. 
One suchphenomenon is the formation of matter-wave 
solitons 0, H M EMMS El El E3 Experi- 
mentally, dark solitons [3, LLfJ and bright gap solitons 
|lfjj were observed in BECs with repulsive interactions, 
whereas bright solitons [l3lll^ were demonstrated in sys- 
tems with attractive interactions. These experimental 
results are augmented by extensive theoretical work in- 
cluding predictions on bright 0, B and dark E3 matter- 
wave solitons, lattice solitons |l2|, and soliton trains 
El- The vast majority of previous theoretical efforts on 
matter-wave solitons have utilized the zero-temperature 
GP mean-field theory. However, in a realistic system, el- 
ementary excitations arising from thermal and/or quan- 
tum fluctuations are always present E3j an d the BEC 
dynamics may be considerably affected by the motion 
of excited atoms around it (thermal cloud), and by its 
dynamical depletion [Isj . giving rise to new nonlinear 
matter-wave phenomena. 

Here we analyze these aspects of BEC soliton dynamics 
by using the time-de pen dent Hartree-Fock-Bogoliubov 
(TDHFB) theory El THE! ■ Soliton dynamics is an- 
alyzed by first calculating the finite-temperature ground 
state of the attractively interacting gas within an har- 
monic trap E3l- The harmonic confinement is then sud- 
denly turned off, and the partially condensed Bose gas 
starts to dynamically evolve. Within the TDHFB model, 



we find a characteristic pattern of evolution of the sys- 
tem whereby pairs of atoms are collisionally excited from 
the BEC into the thermal cloud causing the initial den- 
sity to eventually split into two solitonic structures with 
opposite momentum. Both solitons constitute a mixture 
of condensed and noncondensed particles. We emphasize 
that the observed composite waves are a truly novel type 
of matter-wave solitons, where localization is attained 
not only in spatial density, but also in spatial correla- 
tions. This type of solitons are reminiscent of composite 
incoherent optical solitons 00E3, thus highlighting 
an analogy between incoherent light behavior in nonlin- 
ear media and BECs at finite-temperatures. 

Starting with a near-unity condensate fraction at very 
low temperatures, the GP dynamics reproduces, under 
proper conditions, the well-known zero-temperature BEC 
solitons, thus demonstrating the condensate's mechani- 
cal stability. However, the evolution of the same initial 
nearly pure BEC with TDHFB clearly illustrates BEC 
depiction through pairing, causing these coherent soli- 
tons to disintegrate in a characteristic fashion into in- 
coherent solitary matter waves. In all cases, when both 
the trap and the interparticlc interactions are turned off 
simultaneously, we observe fast matter-wave dispersion. 

We consider a system of N interacting bosons placed 
in a quasi one-dimensional (Q1D) cigar-shaped harmonic 
potential V ex t(x,y, z) = (ui x x 2 + u±y 2 + uj±z 2 )/2, where 
u>j_ 3> lo x denote the transverse and the longitudinal fre- 
quencies of the trap, respectively. The interparticle in- 
teraction is approximated by the Q1D contact potential 
V(xi - x 2 ) = giD$(xi - X2 )) where g 1D = -2h 2 /ma 1D , 
am ~ — a\/azD is the effective ID scattering length 
HI HI E3- TO is the particle mass, a± = yfh/muj±_ is 
the size of the lowest transverse mode, while ol^d is the 
3D scattering length. At finite temperatures, the equi- 
librium state of the system can be described by the HFB 
theory E3 : 
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Here, H sr> = —-^-^j + \mJix 1 and u is the chem- 
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ical potential. The superscript s denotes the static 
HFB calculation, e.g. the static order parameter is 
denoted by $ s (x) and n s c (x) = |$ s (x)| 2 denotes the 
static HFB condensate density |2^| . The normal density 
is n s (x) = Y*j\ uS 3 {x)\ 2N j + \v s 3 {x)\ 2 {N 3 +1) , whereas 
rh s (x) = —J2j u j( x ) v j*( x )(2Nj + 1) is the anomalous 
density |l7l | . The population of excited states at tempera- 
ture T follows the Bose distribution N = (e E ^ kT - 1) -1 . 
In Eq. © we denote C s (x) = H sp + 2gi D [n 6 c (x) + 
n s (x)} - \i and M s {x) = -giD[^ s2 (x) + m s (x)}. Two 
variants of the HFB formalism are the Hartree-Fock 
(HF) approach, altogether neglecting anomalous terms, 
and the HFB-Popov approximation wherein the non- 
condensate anomalous density is explicitly dropped, i.e. 
M s (x) w -g 1D >S> s2 (x) [13. Unlike the static HFB, the 
HF and HFB-Popov approximation do not have an un- 
physical gap in their excitation spectra 0| ■ 

We use the solution of Eqs. and J5J) as the ini- 
tial conditions to study dynamics without confinement 
within the TDHFB approximation in the modal form [2lJ 
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where C(x,t) = H sp + 2gii}[n c (x,t) + n(x,t)], 
M(x,t) = -g 1D {$ 2 (x,t) + rh(x,t)], n(x,t) = 
^Zj\uj\ 2 Nj + \vj\ 2 (Nj + \) is the normal density, and 
rh(x,t) = —^ZjUjV*(2Nj + 1) is the anomalous den- 
sity; at t = 0, all dynamical quantities are identi- 
cal to their static counterparts, e.g., $(x,0) = $ s (x) 
etc. The TDHFB model is usually given in the form 
of coupled equations for the condensate order param- 
eter and the single particle density matrix (e.g. see 
Ref. H3). A tedious but straightforward calculation 
[2^ | shows its full equivalence to the modal form Q). As 
expected, when the time-dependence of the BEC and 
qusiparticle functions is §>(x,t) = § s (x) exp(— ifit/h), 
Uj(x,t) = Uj (x) exp(— i(£"j + (i)t/H), and Vj(x,t) = 

(x) exp(— i(Ej — fx)t/H), the equations of motion (J3J) 
and (J2J reduce to the time-independent HFB equations 
and J2J. 

In what follows we present numerical results based on 
the described formalism, demonstrating the effect of ther- 
mal particles and condensate depiction on matter-wave 
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FIG. 1: (color online) (a) Total density (solid line), and con- 
densate density (circles) at t = 0; total (dot-dashed line) 
and condensate density (dashed line) at tu)± = 180. (b) The 
dipole-shaped first excited state, u\{x) (solid line) and v{(x) 
(dashed line), (c) Evolution of the density with interactions 
present, and (d) without interactions, (e) The complex degree 
of coherence fx(x, x' , t) of a partially-coherent matter-wave at 
t = (dot-dashed line), and tw± = 180 [Re fi(x,x',t) solid 
line, and Im n(x,x',t) dashed line]; x' is placed at the posi- 
tion of the left peak, (f) The condensate fraction (solid line) 
and the total population within the first 60 excited states 
(dashed line) during evolution. 



soliton dynamics. The parameters of the calculation are 
chosen to resemble the experimental parameters of Rcf. 
[3. We consider N = 2.2 10 4 7 Li atoms in a harmonic 
trap with uj± = 4907 Hz (a± = y/h/muj± « 1.35 //m). 
The 3D scattering length cl^d = —3.1 10~ n m corre- 
sponds to a nonlinear parameter of iV|a3rj| « 0.68 urn, 
and is tunable by the Feshbach resonance technique [l3j . 
We present two calculations for different T and lu x to 
show the influence of excitations on the soliton dynam- 
ics. 

First we consider the gas at higher temperature, where 
the condensate fraction is approximately 25%. The trap- 
ping frequency is u)± = 141 Hz (a x = ^Jhjmuj x « 
8.0 /im), while ksT/hw^ = 60. Fig. da) illustrates the 
total density n 8 c (x) + nf(x) of the stationary HFB-Popov 
calculation. The total density profile is double humped, 
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which is a clear signature of the significant population of 
the first excited state (approximately w 10% of atoms), 
which has a dipole-like shape shown in Fig. IHb). We 
have numerically checked the stability [3(| of the solu- 
tion with respect to small perturbations; the stability 
is underpinned by the use of parameters resembling ex- 
periment When the trap is turned off, the system 
is suddenly taken out of equilibrium, and consequently 
starts to evolve; we simulate the dynamics with the full 
TDHFB model. In the spirit of Ref. 0|, we compare 
the x-unconfined dynamics of the system in the pres- 
ence of intcrparticle interactions [Fig. dc)] to its time 
evolution when both the confinement in x and the in- 
teractions are turned off [Fig. IHd)]. In the absence of 
interactions, we clearly observe fast matter-wave disper- 
sion. In contrast, when interactions are present, the two 
humps begin to separate, because the trapping poten- 
tial which provided a balance to the kinetic energy term 
is no longer present [Fig. dc)]. During evolution, the 
condensate is slowly depleted [Fig. IHf)]. Consequently, 
at one point during the evolution the two humps split, 
each forming a solitonic structure with opposite momen- 
tum. Each of these solitonic structures contains both 
Bose-condensed atoms, and a significant portion of the 
non-condensed particles, i.e., they are partially coherent 
matter- wave solitons. The uniqueness of such random- 
phase structures is ellucidated by their complex degree of 
coherence, p,(x,x',t) = p(x, x' , t)/y/ p(x, x, t)p(x', x',t), 
where p(x,x',t) = $(x', t)*$(x, t) + (&(x',t)$>(x,t)), 
plotted at different times in Fig. [He). We observe 
that spatial correlation is localized in space at t = 0. 
Correlations change as the two humps split. While re- 
taining spatial localization, the phases at the two sep- 
arated peaks are well-correlated, and out of phase, re- 
sembling the behavior of out-of-phase adjacent solitons 
from Ref. [l4|. We emphasize that for zero-temperature 
GPE solitons, the pair correlation function factorizes as 
p(x,x') = &*(x)$(x'), which yields p(x,x') = 1, cor- 
responding to coherent matter-waves. Our incoherent 
matter-wave solitons are thus rather special in that they 
correspond to localization of entropy and spatial correla- 
tion, as well as to localization of density. 

While the temperature fc^T in the previous example is 
higher than the transverse level spacing fku i whereas a 
'true' ID geometry calls for ksT < Hu)± [23, the use 
of a Q1D formalism is still justified because the first 
u> x /u)± ~ 35 states are essentially ID (they are in the 
lowest state of the transverse Hamiltonian) . Further- 
more, only condensed atoms, and atoms from lower ex- 
cited states determine the outcome of the motion. There- 
fore, a proper inclusion of the transverse dimension in the 
calculation would lead to some rescaling of the parame- 
ters, but would not influence the dynamics observed in 
our quasi-lD calculation. Moreover, the simulations as 
well as the experiment of [l3T |. are all in the weak inter- 
action regime N\aio\/ax ~ 10 8 3> 1 H3i thus justifying 




-20 tm 20 40 %, 




20 w'a 


20 


40 


.(d) ;\ 















-20 X&. 20 40 



FIG. 2: (a) The evolution of the system with GPE, and (b) 
with TDHFB. (c) The condensate fraction (solid line) and 
the fraction of non-condensed atoms during TDHFB evolu- 
tion, (d) The complex degree of coherence n(x,x',t) at t — 
(dot-dashed line) , and after the condensate is considerably de- 
pleted at tu>± = 73 [Re n(x, x' , t) solid line, and Im p(x, x' , t) 
dashed line]; x' is placed at the position of the left peak. 



the use of a mean-field approach. 

Next, we consider a gas prepared at near-zero tempera- 
ture, where the condensate fraction is 99%. The trapping 
frequency is uj± = 439 Hz (a x = y/h/muj x w 4.51 /im), 
while kBT/hui± = 5. As previously, the trap is turned 
off and the dynamics is calculated using both the GP 
equation (Fig. El) and the TDHFB formalism (Fig. OJd). 
The mechanical stability of the evolving BEC is demon- 
strated by propagating the time-dependent GP equation, 
showing small oscillations of the condensate width, rather 
than a mechanical collapse to a point [Fig. P[a)]. This 
motion corresponds to a coherent matter wave soliton. 
However, as clearly evident from propagation of the TD- 



HFB equations [Fig. Efb)], allowing BEC depletionjFig. 
Efc)] , the BEC collapses via a pairing instability |3ll l32| 
whereby pairs of atoms are collisionally pulled out of the 
BEC into the thermal cloud, thus gaining a mean-field 
energy which goes into their relative motion. Such pair- 
ing collapse with little or no mechanical shrin king is in- 
deed observed in 3D collapse experiments j^S |3| ■ Our 
results show that in ID it results in two separate solitonic 
structures of opposite momentum, containing both a con- 
densed part and a significant thermal population. Wo 
note that similar structures were observed in stochastic 
simulations of molecular BEC dissociation in ID geome- 
try 1^3, indicating that incoherent matter- wave solitons 
can also be produced in this system. Figure H^d) shows 
that the correlations of each solitonic structure are local- 
ized, and the two structures are partially correlated. 

Before closing, we note that within the time-dependent 
Hartree-Fock (TDHF) approximation where con- 
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dcnsatc depletion is not accounted for during the dynam- 
ics, we find solutions with double-humped density profile 
that begin to split as in TDHFB model. However, as the 
condensate is not depleted within TDHF, the two peaks 
may be pulled back and merge to almost recover the ini- 
tial density profile. Such motion is also characteristic of 
composite incoherent solitons in optics. 

In conclusion, we have used the time-dependent 
Hartree-Fock-Bogoliubov theory to analyze the influence 
of the thermal cloud and condensate depletion onto the 
dynamics of BEC solitons. We find that condensate de- 
pletion induced by pairing, and the presence of a thermal 
cloud cause the particle density to split into two solitonic 
structures, each being a mixture of condensed and non- 
condensed particles. The predicted incoherent matter- 
wave structures represent novel correlation solitons which 
resemble localized second-sound entropy waves. Such 
random-phase matter-wave solitons correspond to inco- 
herent solitons in nonlinear optics j^, Ej, El- which 
points at the analogy between partially condensed Bosc 
gases and nonlinear partially coherent optical waves. 
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